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We investigate the issue of speed-up and the necessity of entanglement in Grover's quan- 
tum search algorithm. We find that in a pure state implementation of Grover's algorithm 
entanglement is present even though the initial and target states are product states. In 
pseudo-pure state implementations, the separability of the states involved defines an en- 
tanglement boundary in terms of a bound on the purity parameter. Using this bound 
we investigate the necessity of entanglement in quantum searching for these pseudo-pure 
state implementations. If every active molecule involved in the ensemble is 'charged for' 
then in existing machines speed-up without entanglement is not possible. 
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1. Introduction 

It was Feynman's insight to realize that quantum systems cannot be efficiently simulated 
on conventional classical computers Subsequently, Deustch suggested that if one could 
build computers exploiting the principles of quantum theory one might be able to speed up the 
computation process compared to classical approaches ■ Quantum computers aim to make 
use of quantum interference and entanglement between different parts of a bulk quantum 
system to give such an essential difference ||. 

In recent years an important question has been: What will it take for quantum computers 
to surpass their conventional classical counterparts in speed and efficiency? The answer 
consists of algorithms and their efficient implementation. The first promising algorithms for 
quantum computers were discovered by Deutsch and Jozsa for function testing by Shor 
for factoring Q and by Grover for searching |J. While the efficiency of these algorithms 
is today well established, the conditions for achieving quantum efficiencies have been the 
subject of recent controversy |7). Indeed, experimental implementations of Grover's algorithm 
on pseudo-pure state machines using liquid-state Nuclear Magnetic Resonance (NMR) || ^] 
have been claimed to already achieve quantum efficiencies JTc| , [Tl| , in spite of their 

apparent inability to have produced entanglement to date Similarly, the Deutsch- Jozsa 
algorithm has been implemented on NMR quantum computers [^5[ 

The aim of this paper is to argue that the original version of Grover's algorithm H 
on multiple qubits necessarily involves quantum entanglement, even though the initial and 
target states are product states. Further, by counting each active molecule as contributing 
to the computational resources for pseudo-pure state machines, we show in a non-asymptotic 
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analysis that not only is entanglement necessary to achieve a speed-up in quantum searching, 
but it must be present throughout the computation. If one uses a different resource counting 
method there may be speed-up without entanglement. Thus, we have resolved the role played 
by entanglement in quantum searching on one whole class of quantum computers. This non- 
asymptotic result unequivocally proves that a class of liquid-state NMR machines jlO], [LlJ [l^, 
|f3| cannot be used for performing faster-than-classical quantum computation. 

2. Efficiency in Quantum Computation 

Typically, efficiency is quantified by relative 'speed' or how the number of steps needed to 
complete an algorithm scales with the size of the 'input' the algorithm is fed. Two ubiquitous 
'exponential' problems are searching and factoring: All known algorithms for solving them on 
conventional computers scale roughly exponentially with input size (e.g., the length of the list 
to be searched or size of the number to be factored). Discoveries of 'fast' quantum algorithms 
j| |[ D set new bounds on computational goals and standards. 

In order to determine the efficiency of an algorithm on a quantum computer, the conven- 
tional measure of 'speed' must first be re-evaluated. Clearly, for the scaling behavior to be a 
sensible measure of efficiency that may be used to compare the performance of very different 
kinds of computers, there must be no 'hidden costs' that grow in an unreasonable manner 
(i.e., faster than the scaling itself). For example, any increase in the size of the computer 
itself, or number of resources it utilizes relative to the input, should not exceed the scaling of 
the number of steps. 

Pseudo-pure state quantum computers are welcome candidates for such re- valuations. In- 
deed, such machines have been implemented using liquid-state NMR and have been proposed 
for a wide variety of quantum algorithms (|, ||, [l5|, |l(| [l7| and their efficiency has already 
been determined for asymptotically large systems (i.e., many qubits) for Shor's factoring al- 
gorithm. It was found that, in this asymptotic limit, an absence of entanglement would lead 
to an exponential decrease in the probability for obtaining the correct answer [ tLq |. Thus, if 
Shor's algorithm were to be implemented on such machines, an exponentially large number of 
resources would be required to boost this low probability. Because Shor's algorithm provides 
an exponential speed-up, it may not be so surprising that the weirdest features of quantum 
mechanics, i.e., entanglement, are required for its implementation. 

An important caveat of this analysis stems from the extreme unlikelihood of constructing a 
sufficiently large machine that could be subjected to such asymptotic analysis. For example, 
it was observed that for liquid-state NMR implementations the signal scales as n/2 n with 
increasing numbers n of qubits (i.e., active spins) Jl9| . This inherent scaling problem suggests 
that regardless of entanglement or speed-up there seems little hope of ever reaching the 
asymptotic regime. 

The situation is very different for Grover's search algorithm. First, Grover's quantum 
algorithm provides a much more modest quadratic (as opposed to exponential) speed-up 
over any search on a conventional computer; thus one might expect it to be more robust 
with respect to a loss of entanglement. Secondly, there have now been several experiments 
demonstrating this algorithm on small (few-qubit) NMR machines and claiming quantum 
efficiencies [jl(J [fl], |l2|, 13 . For these experiments, the asymptotic signal scaling mentioned 



above is not relevant, and the entanglement question thus becomes the primary consideration 
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in assessing speed-up. 

In fact, speed-up in Grover's algorithm can be evaluated quite naturally in terms of query 
complexity. The query complexity formulation yields a non-asymptotic result which may 
be applied to any size problem. We argue that, in pure-state as well as pseudo-pure state 
implementations, if only separable (i.e., unentangled) states are accessed then the speed-up 
predicted by Grover's algorithm fails to materialize. There exists a peculiar exception (though 
one not accessible to existing pseudo-pure state machines) for a search space of size four where 
entanglement is not necessary. Since generalizations of this interesting exception are as yet 
unknown, the heuristic claim that entanglement is necessary for a scalable quantum compu- 
tation implementing Grover's algorithm still holds. Our conclusion is also consistent with 
the analysis of Grover's algorithm, implemented on a device exploiting only superposition, 
but at the cost of scalability due to exponentially growing resources Q . It should be noted 
that in addition to the necessary condition we analyze the sufficient condition for speed-up of 
quantum searching on pseudo-pure state machines. 

The outline of the paper is as follows. In section 3, we begin with a straight-forward 
demonstration that the pure-state version of Grover's algorithm involves entanglement during 
the quantum searching operation. In section 4, we analyze the quantum search that can be 
implemented on pseudo-pure states and derive a bound on the purity parameter for the 
pseudo-pure states to be separable. In section 5, to determine the necessity of entanglement 
in the pseudo-pure state version, two complementary criteria are derived: one for the presence 
of entanglement (as a function of the number of qubits) and one for the query complexity of 
the algorithm, yielding a measure of speed-up relative to the classical algorithm. We impose 
these two criteria in the few-qubit regime, one at a time, and obtain a one-to-one relation 
between entanglement and speed-up, showing both that there can be no speed-up without 
entanglement, and, conversely, that in case a speed-up is in fact achieved, entanglement must 
be present throughout the computation. In section 6, we discuss how pseudo-pure states are 
not good candidates for mimicking pure-state dynamics. Finally, we discuss the implications 
of our results. 

3. Quantum searching with pure states and entanglement 

According to the standard formulation of the search problem || , we are given an unknown 
binary function f(x), which returns 1 for a unique 'target' value x = y and otherwise, over 
the domain x = 0, 1, 2, . . . , N — 1, with N — 2". Our goal is to find y such that f(y) = 1. 
In Grover's algorithm, the N inputs are mapped onto the states of n quantum bits (qubits) 
such as spin-i particles. The quantum problem thus becomes one of maximizing the overlap 
between the state of these n qubits and target state \y). This is equivalent to maximizing the 
probability of obtaining the desired state upon measurement. The initial state of these qubits 
is taken to be an equal superposition of all possible bit strings, i.e., 



The Grover operator defined as G = —Jo H® n l y H® n is used repeatedly in the algorithm, 
where Iq — 1 — 2|5 , o)(^ / o|i ly = 1 — 2|y)(j/| with \y) being the target (ideally the final) state 
and H is the Hadamard transformation. Thus, the Grover operator corresponds to a small 
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rotation in the two-dimensional subspace spanned by the initial and target states. Each such 
rotation requires a single evaluation of /. Thus, unlike a classical search, the quantum search 
monotonically rotates the state towards the target. 

Let us start by considering the pure-state version of Grover's algorithm. After k iterations 
of the Grover operator the combined n qubit state (|l|) evolves to j2lj 

|* fe ) = -^=^|x)+sm«, (2) 



where 9k = (2k + 1)6*0 and 6q satisfies sin 6*0 = 1/\N. The search is complete when Ok — ft/2 
which takes 0(-J~N) iterations of the Grover operator and hence this many evaluations of the 
function /. 

We will now show that although the initial and target states are product states the interme- 
diate states through which system evolves are entangled. Since these states are superpositions 
of product states they are expected to be entangled. But how much entanglement is there 
in these intermediate states? This is difficult to quantify, as we do not have a proper mea- 
sure of entanglement for quantum systems consisting of an arbitrary number of subsystems. 
However, we can consider the system as being bipartite, with one subsystem consisting of a 
single qubit and the second subsystem all the rest. In this way, we will be able to quantify 
the bipartite entanglement by calculating the reduced density matrix of any single qubit. 
We then ask how close this reduced state is to a maximally entangled qubit (using say the 
Hilbert-Schmidt norm criterion). 

The reduced density matrix of the fth qubit is 

ptQt) = tT 1>2t ..,e-l,i + l,.. n (\^k}{^k\) 

= alH\Q)(m + bl\jd(jA + ^m)(3i\ + \u)(jA + \u)QA) , (3) 



where ak = \JN/(N — 1) cosflfc, bk — sin 9k — cos 8k /^/N — 1 and the single bit jt = 1 — ji, 
(je — 0, 1). Without loss of generality we take jg = 1 and the density matrix pt(k) can be 
expressed in standard form, i.e., 

Pl (k) = \[i + m •<?] = [!- *(*)]§ + s(k)p, (4) 

where s(k) = tr \pi (k)a], s(k) ■ s(k) — s(k) 2 < 1 and P is a pure state projector. The 
components of the Bloch vector s(k) after k iterations are 

s x (k) = cos 2 9k H sin 9.9, 

y ' N-l y/W^l 

Sy (k) = 

s z (k) = — -j- cos 2 9 k - sin 2 9 k ■ (5) 

The bipartite entanglement in the pure state may be characterized by calculating the von 
Neumann entropy of this reduced state. Using the expansion formula 

log(^l + sP) = log(^) + P log(l + ^) , (6) 
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which holds for any < s < 1, the von Neumann entropy may be calculated to be given by 
S\p t (k)] = -tr [ Pe (k) log p e (k)} 

= l- 1 —^\o g [l-s(k)]- 1 -±^\og[l + s(k)}. (7) 

The right-hand-side of this expression is independent of the choice of the remaining qubit 
I. Therefore, (7) holds for any one qubit versus (n — 1) qubit partioning. It shows that the 
reduced density matrix of the single qubit does not arise from a maximally entangled state of 
n qubits, as the von Neumann entropy is not exactly unity. Since the reduced state of Eq. (|^) 
is not pure the full state must be entangled. To see how impure the state in Eq. (^) is one 
may calculate the linear entropy L{p) of it which is given by 

L[p t {k)] = tr \p t (k) - Mk) 2 } = 1 ~ * (fc)2 • (8) 

If the linear entropy is zero the state is pure and as it approaches i the state approaches 
a completely random mixture. In the quantum search algorithm the parameter s(k) can 
never be zero because that would mean that cos 9k and sin 9k are simultaneously zero, which 
cannot be satisfied. So although the reduced density matrix of the qubit may lie close to the 
completely mixed state it can never become the identity one. 

We now calculate the Hilbert-Schimdt norm of the difference of two density matrices: 
namely, the completely mixed one and the reduced state given by Eq. (|J). This will yield a 
measure of closeness to the completely mixed state. This Hilbert-Schmidt distance for fcth 
iteration during quantum search algorithm is given by 

d(kf = \\~-p t (k)\\ 2 HS =tr[~-p t (k)} 2 

= l-L[Mk)] = ^. (9) 

The distance d(k) provides an idea of how the reduced state of an individual qubit behaves 
during the fcth iteration. It shows that the reduced density matrix of the qubit differs from 
a completely random mixture by an order of 0(s(fc)). One can see from Eq. (§) and (§) that 
for 9q = sin (1/ViV) and for 9k — 7r/2 the reduced density matrix of any remaining qubit 
is pure, implying that the whole state must have been non-entangled. Thus, we see that 
although the initial and target states are separable, the intermediate states through which 
the system evolves are always entangled. 

4. Quantum searching with pseudo-pure state and entanglement 

In the following we discuss the issue of entanglement in quantum searching with pseudo- 
pure states. These states naturally arise in liquid-state NMR machines. Here one faces the 
difficulty of accessing a pure state because the system is in thermal equilibrium at room 
temperature. Instead, one implements Grover's algorithm on a near random ensemble of 
molecular spins in a liquid, with a small preference for the spins to point along an external 
magnetic field; the size of this preference is quantified by the purity parameter e [typically 
O(10~ 5 )] H ||, [l^, [H], |l^, ^3|. For a sufficiently low spin polarization (corresponding to a 
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sufficiently low purity parameter), the system can be well-approximated by a pseudo-pure 
state representation described by 

p =izfi JV + e |*)(*|, (io) 

where 1 at is the identity matrix of dimension N. (In the high e regime, an exact description 
can be obtained by reverting to the full Boltzmann distribution.) Since, for experimentally 
feasible values of e (now and in the foreseeable future), the pseudo-pure state description 
suffices, we henceforth restrict our discussion to them. 

Let us now consider quantum searching on pseudo-pure quantum states. After k iterations 
of the Grover search operator, one obtains the state 

Pk = ^ 1 lN + e\H> k )(V k \ , (11) 

where \^f k ) is given by Eq. (||) 

Let us first look at the diagonal form of the reduced state of any remaining qubit (denoted 
by the index £) in the pure state given in Eq. (|1). This has positive eigenvalues Ai(fc), A2(fc) 
which are independent of £ but depend on the iteration index k. These eigenvalues sum to 
one and their product, after some calculation, is given by 



Ai(fc)A 2 (fe) = ^ sin 2 (2fc0 o ) cos 2 6 k . (12) 

This allows us to decompose the full state at step k into a Schmidt basis 



|* fc ) = VMQWM - Vxm\e')\g) , (13) 

where {|e), \g)} describes an orthonormal basis for the £th qubit and {|e'), \g')} is a pair of 
Hilbert space vectors for the other (n — 1) qubits. This expression will help us to derive a 
bound on the purity parameter. 

In order to study the entanglement present in the pseudo-pure state given by Eq. ( |ll| ) 
let us project it onto the 4-dimensional subspace spanned by the set {\g')\g), \g')\e), |e')|<?), 
|e')|e)}. The resulting 4-dimensional density matrix is 

^-4T7(^i)(V l4 + ^^')- ^ 

The boundary between separability and entanglement for such 4-dimcnsional states may 



be specified in terms of the so-called fidelity F. The states (14) are separable when F 



(^~\pi^\^~) - \ and are entangled otherwise |22j], where |* _ ) cx \g')\e) - |e')|g). Since 
projection onto a subspace cannot create entanglement, it follows from this condition that if 
the original unprojected state (|ll]) is separable, then we must have 

e < e fc = 1 . (15) 

l + AVAi(fc)A 2 (fc) 

Similarly, the original state (at stage k of the search) is entangled whenever e > e k - Similar 
bounds have been earlier derived when the pure state part was a maximally entangled state 
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p4[ . However, the present bound quantifies the separability region of pseudo-pure states for 
each iteration k during Grover's algorithm. Thus one can tell at each stage of computation 
whether the state is entangled for a given value of purity parameter. 

5. Speed-up and Entanglement for Pseudo-pure states 

Let us now return to a careful re-evaluation of what it means for a pseudo-pure state 
implementation of Grover's algorithm to demonstrate a speed-up compared to a classical 
search on a conventional computer. First, we note that the important common element 
in either a quantum or a classical search algorithm is the repeated need of evaluating the 
function /. Thus, we shall restrict our comparison to its number of evaluations. This is 
formally known as studying the query complexity of the problem. Second, since Grover's 
algorithm is probabilistic, there is no upper-bound to the worst-case number of evaluations. 
Instead, we study the expected number of evaluations for finding the target. To be fair then, 
we must compare this with the expected number of function-evaluations on a conventional 
computer. If we exclude the use of an exponentially growing auxiliary memory to store failed 
trials then this classical query complexity is 

Ai_ . (" + »«"-'> PA 

for finding a single object out of N entries. This value may be achieved by systematically 
stepping through the N possible locations for the object and evaluating / to see if the object 
is there. If it is not found in this way by step N — 1 then we would know it is at the final 
location N. (The specific value for classical query complexity quoted in Ref. jl3| for N = 8 
can be obtained from our general result 7V c iass above). 

How does this compare to the expected number of function evaluations for quantum search 
on pseudo-pure state implementations? The probability of finding the target state after k 
iterations is 

/ , , x fl + e(7Vsin 2 6> fe - 1)1 , . 

p(k) = (y\p k \y) = 1 v N - < 1. (17) 

This probability must be amplified statistically through repetitions or parallelism. The ex- 
pected number of repetitions required to identify the target is just the reciprocal of this 
probability. Each such repetition involves k function-evaluations for each run of the algo- 
rithm, plus one to test the result. In order to give Grover's algorithm its maximal advantage, 
we shall optimize the speed (rather than the probability). Thus, the optimal expected number 
of function evaluations is just 

_ min fc (fc + l) _ (fcopt + 1) n „s 
A/ pscudo - p{k) - , U»j 

where k op t is the optimal number of iterations of the Grover operator. Thus, a pseudo-pure 
state quantum computer can search faster than a conventional computer, provided 

-^pseudo ^ -^class • (1^) 



In the above analysis we have utilized the standard projective measurement which is a 
desirable operation in any quantum information processing unit. However, the liquid-state 
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NMR machines use weak measurements and the success probability obtained in a strong 
measurement is always greater than the success probability for a weak measurement. So in 
our analysis we give extra benefit to the experimental situation in tolerating the error. In 
any case, if we take weak measurement the situation for NMR machines will be no better (as 
that will increase the number A/p Se udo)- 

We have now developed the tools to effectively answer the question: Can Grover's algo- 
rithm yield a speed-up on any existing pseudo-pure state implementations [fLO] , [j~lj [l2[ fL3| ? 
Since all liquid-state NMR experiments performed so far have generated only separable en- 
sembles fl4|| , with purity parameters that fit well our pseudo-pure state treatment above, 
let us see what effect separability has on efficiency. From Eq. (|T^) separability throughout 
the search implies e < min^p e^. Thus, separability places an upper bound on the purity 
parameter (see values for this bound in Ref. [Q) and hence a lower bound, A/^^^ , on the 
quantum query complexity. This lower bound is given in Table I for the optimized Grover 
algorithm for search spaces up to n = 8 qubits. (The trend we find continues for arbitrarily 
large numbers of qubits.) The surprising observation is that, for more than 4 qubits, the 
optimal quantum strategy is not to use the liquid-state NMR machine at all and simply guess 
the answer (equivalent to sampling the initial random state). Clearly, this is never as good 
as a systematic conventional search. 



Table 1 



n 


N 


k pt 


.r(min) 

pscudo 


•A/class 


1 


2 





2 


1 


2 


4 


1 


2 


2.25 


3 


8 


1 


5.48 


4.38 


4 


16 


2 


12.89 


8.44 


5 


32 





32 


16.47 


G 


64 





64 


32.48 


7 


128 





128 


64.49 


8 


256 





256 


128.50 



Curiously, we note that our table shows a speed-up even without entanglement for the 
n = 2 qubit implementation. That this algorithm requires no entanglement in this case, has 
also been noted for the pure-state implementation of Grover's algorithm |20(] (with a similar 
result for the two qubit implementation of the Deutsch-Jozsa algorithm fl23| ). This occurs 
because the separable target state can be reached by a single application of the Grover op- 
erator from the separable initial state; hence the evolution has passed through no entangled 
states. Without entanglement in the two-qubit case there is also no penalty to the efficiency 
when enforcing separability through a reduced purity parameter. Despite this curiosity, this 
speed-up is only possible for unphysically large purity parameters, e > 23/27 ~ 0.852, where 
the pseudo-pure description is no longer valid for liquid-state implementations. This obser- 
vation and the above analysis show that entanglement is necessary for obtaining a speed-up 
in Grover's algorithm on a liquid-state NMR machine relative to a classical computer. Let 
us emphasize that for pure-state machines (e.g., ion-traps) the two-qubit implementation of 
Grover's algorithm is realistic and will yield a speed-up of 2.25 (omitting the now superfluous 
final testing function-evaluation). 

As we have seen, the presence of some entanglement is essential to obtain a speed-up, 
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except in the two-qubit case. But how much of it? Let us impose the speed-up condition ( |l9| ) 
to obtain a lower bound e spC cd-up(&opt) on the purity parameter within the pseudo-pure state 
approximation. After iteration k, condition (15) implies entanglement is present whenever 
e sp eed-up(fcopt) > £fc- We studied this relation numerically for 2 < n < 20 qubits with 
< k < fc pt < \tt/40o — 1/21, an d found that entanglement was present after every iteration 
except possibly the last step of the algorithm when 0& t > it/2. Thus, we may draw the much 
stronger conclusion that, for more than two qubits, some degree of entanglement is necessary 
during the entire computation in order to obtain any speed-up for Grover's algorithm on a 
liquid-state NMR machine. 

To better understand the repercussions of this result, let us directly address the main 
objections raised about the role of entanglement in performing real quantum computation in 
this system. The most dismissive argument has been that quantum computational efficiency 
derives only from the unitary evolution of quantum states, but is independent of the type of 
states being used (Laflamme in Ref. jjj). Now, for the sake of clarity, let us adopt the defi- 
nition of quantum computing suggested by the same proponent in the same context, namely: 
Quantum computing consists of efficiently evolving from initial to final density matrices by 
unitary operations. Since it is not enough to reach the final state probabilistically, we have 
calculated the query complexity of the problem in detail (assuming perfect unitary evolu- 
tion). Thus our analysis clearly proves that unitarity alone is not enough to achieve quantum 
efficiencies. 

A different objection correctly points out that liquid-state NMR machines are not exactly 
described by the pseudo-pure state formalism, and hence the bounds we derive above for 
entanglement and separability may not be applicable. Indeed, corrections to this bound are 
necessary in the few-qubit regime, where the bounds for the purity parameter are high. How- 
ever, in most likelihood, these corrections will not alter our results and may even strengthen 
them (by raising the lower bound for entanglement). Further, noting that existing liquid- 
state NMR machines are remarkably far from reaching this range of purity parameters, any 
possibility of speed-up is quite out of the question. 



6. Dynamics of Pseudo-pure versus pure states 

Granted that entanglement is needed to achieve quantum efficiencies, it has been argued 
|L8| that the ability of liquid-state NMR machines to physically follow unitary quantum 
evolution qualifies them as efficient simulators (without speed-up) . The reasoning is based on 
the realization that the observable operators are traceless Thus, the expectation value 
of any traceless observable on a pure quantum state is the same, up to a scale factor, as it 
would be on a pseudo-pure state. This is because for a traceless observable O the average in 
pseudo-pure state is tr(/?6) = e(\l/|0|\l/). However, in the quantum world fluctuations play 
as important a role as expectation values. In general, if an observable is traceless its square 
need not be. In this case the quantum fluctuations of an operator on a pure state are not 
equal to those fluctuations on the corresponding pseudo-pure state. For example, if 6 is a 
traceless operator its root-mean-square fluctuations on pseudo-pure states are determined by 

(Ae)2 scudo = e(A6)^ ure + (1 - e)(^p- + e(e) 2 puxe ) . (20) 
As a simple example, consider = dp/de = P — 1/N, where P = It is a traceless 
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operator, but its square is not. The uncertainty of the operator Q on the pure state is zero, 
whereas on the pseudo-pure state it is 



(A6)^ udo = (l-e)(l-l) ^-e(l-^) 



(21) 



In fact, very rarely will it be the case that all the moments (even for the restricted class of 
traceless operators) are related by simple scale-factors to those produced by a liquid-state 
NMR quantum computation. Thus, despite the correct expectation values being accessible 
for traceless operators, we cannot say that liquid-state NMR quantum computers have good 
dynamics. 



7. Conclusion and discussion 

In conclusion we have shown that the original version of Grover's algorithm implemented 
on qubits necessarily generates quantum entanglement during the computation process. Fur- 
ther, we have shown that a quantum computer in a non-entangled pseudo-pure state requires 
more iterations than a classical computer to perform a virtual database search if one uses our 
method of counting the number of queries. Thus, for any existing liquid-state NMR set-ups 
our results preclude any possibility of speed-up in running Grover's algorithm. This conclu- 
sion is based on exact calculations for pseudo-pure state implementations, which show that 
entanglement is essential throughout the computation. Moreover, even a modest reduction of 
entanglement is tantamount to a total loss of speed-up. Despite the decisive results reached 
for pseudo-pure state implementations of this algorithm, we found that entanglement-free 
speed-up would be possible for the two-qubit case (search space of size four) of Grover's 
algorithm on pure-state machines such as ion-trap quantum computers. Our analysis has 
been with respect to a specific non-asymptotic regime and does not find a way for scalable 
computation without entanglement. Finally, our analysis does not resolve how other NMR 
experiments should to be interpreted such as quantum teleportation without entanglement. 

Note added: After we completed this work, a recent article has even implemented the 
quantum search algorithm using only the wave nature of classical Fourier optics p4| which 
involves no entanglement! However, one should not be surprised into thinking that this is 
a contradiction to our findings presented here. What happens with classical devices imple- 
menting the search algorithm is that the number of resources needed increases prohibitively 
with the size of the input. 
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